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摘要 : 研究 发 现 ， 即 使 在 银河 系 长 期 缓慢 演化 的 过 程 中 ， 其 核 球 中 依然 可 能 存在 垂直 的 金属 丰 
度 梯 度 。 该 结果 反驳 了 早期 研究 中 认为 的 银河 系 核 球 中 长 期 缓慢 演化 的 过 程 会 消除 任何 存在 的 
化 学 丰 度 梯度 的 观点 。 通 过 使 用 一 个 模拟 银河 系 薄 盘 的 多 体 数 值 模拟 模型 得 出 该 研究 结果 ， 该 模 
型 能 自治 地 通过 盘 的 棒 不 稳定 性 和 届 曲 不 稳定 性 形成 一 个 盒 状 /花生 状 的 核 球 ， 并 结合 化 学 -动力 
学 模型 ， 分 析 银 河 系 核 球 中 化 学 成 分 的 演化 。 

在 该 演化 模型 中 设 定 了 一 个 初始 的 径 向 金属 丰 度 梯度 (—0.3 dex/kpc)， 通 过 长 期 缓慢 的 演 
化 过 程 ， 其 核 球 部 分 呈现 出 与 观测 相符 的 垂直 金属 丰 度 梯度 。 对 此 一 种 可 能 的 解释 是 ， 星 系 盘 的 
金属 丰 度 分 布 在 长 期 缓慢 演化 的 过 程 中 经 历 了 “两 步 ”演化 过 程 ， 首 先 盘 星 系 初始 半径 较 大 的 贫 
金属 粒子 在 径 向 混合 的 过 程 中 产生 了 更 大 的 径 向 速度 弥散 ， 然 后 由 于 核 球 的 屈曲 不 稳定 性 演化 ， 
垂 向 速度 弥散 最 终 会 保持 与 径 向 速度 弥散 的 一 定 比值 ( 约 0.8)。 这 个 机 制导 致 贫 金 属 粒子 在 演化 
过 程 中 能 占据 更 大 的 垂 向 范围 ， 从 而 形成 核 球 的 垂直 金属 丰 度 梯度 。 
银河 系 核 球 动力 学 模型 通过 简单 的 化 学 -动力 学 过 程 揭示 的 核 球 金属 丰 度 演化 机 制 ， 应 在 包 
含 更 复杂 的 核 球 演化 机 制 的 模型 中 同样 存在 。 此 外 ， 还 通过 化 学 动力 学 模型 讨论 了 核 球 观测 中 存 
在 的 争议 ， 如 金属 丰 度 分 布 函数 中 产生 多 个 “ 峰 ” 的 起 源 问题 ， 以 及 该 模型 构建 中 存在 的 一 些 的 
不 足 。 
X 键 词 : 银河 系 ， 银 河 系 核 球 ， 银 河 系 演化 ， 银 河 系 丰 度 
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上 一 
E 
li 


近年 来 ， 银 河 系 中 心 核 球 的 观测 和 研究 有 助 于 我 们 到 
中 发 现 的 银河 系 中 心 结构 分 布 和 运动 学 特性 主要 包括 形态 上 不 对 称 的 平行 


41 卷 


解 银河 系 的 形成 和 演化 过 程 。 观 测 
四 边 形 ”、 运 动 学 


上 的 柱状 旋转 ”， 以 及 核 球 恒星 分 布 的 X 型 特征 ” 等。 这些 证 据 揭示 了 银河 系 的 核 球 起 源 


于 星系 盘 结构 本 身 的 长 期 动力 学 演化 。 利 用 数值 模拟 的 银 ; 
河 系 核 球 的 盒 状 /花生 状 结构 从 星系 盘 中 天 


能 较 好 地 匹配 银河 系 核 球 的 观测 结 
模拟 银河 系 核 球 更 细致 的 结构 时 ， 


发 现 银河 系 


"zi 


中 心 并 不 包含 一 个 显著 的 经 典 核 


bulge)。 这 意味 着 ， 银 河 系 更 接近 于 一 个 纯粹 由 星系 盘 通 过 自身 的 长 期 缓慢 演 


Anm 


不 过 也 存在 一 些 和 否认 银河 系 长 


盘 平 面 的 金属 丰 度 梯度 ， 即 随 着 银 纬 增加 核 球 


动力 学 演化 总 是 涉及 到 棒状 不 稳定 


在 银河 系 的 动力 学 长 期 缓慢 演化 的 


的 构建 来 理解 银河 系 核 球 中 的 垂 向 金属 3 


期 演化 图 景 


的 观点 ， 如 在 银河 系 核 球 中 观测 到 的 垂直 于 
区 平均 金属 丰 度 下 降 的 现象 。 由 于 银河 系 长 期 


性 和 屈曲 不 


Bekki 和 Tsujimoto AI, ZEA Av 


系 模 型 ， 许 多 工作 能 够 还 原 出 银 
E 成 的 过 程 ， 并 且 其 关于 动力 学 和 形态 学 的 结果 也 
ARP". 2010 ££, Shen 等 人 "提出 银河 系 核 球 模型 ， 在 


球 (classical 
化 形成 的 星 


稳定 性 演化 的 强烈 动力 学 混合 过 程 ， 不 少 天 
文学 家 认为 这 种 演化 不 太 可 能 形成 一 个 有 序 的 化 学 丰 度 梯度 ””。 为 了 解决 这 个 问题 ， 需 要 


化 模型 中 ， 当 给 予 初始 盘 核 球 3 


图 景 中 加 入 化 学 丰 度 的 分 布 情况 ， 通 过 化 学 -动力 学 模型 
FE 度 分 布 的 形成 和 演化 过 程 。 


径 范围 ( 约 


2 kpc) 内 的 径 向 和 垂 向 金属 丰 度 梯度 的 分 布 后 ， 确 实 出 现 了 演化 后 形成 的 核 球 中 垂 向 金属 梯 
度 较 平 的 情况 。 随 后 Martinez-Valpuesta 和 Gerhard 守 指出 ， 要 实现 符合 观测 
重 向 金属 丰 度 梯度 的 分 布 ， 初 始 星 系 盘 所 有 涉及 到 棒 / 核 球形 成 的 范围 内 ( 即 棒 共 转 半径 以 


内 ， 约 4.5 kpc) 都 要 存在 一 定 的 金属 丰 度 梯度 。 运 用 这 个 方法 ，1 


中 的 银河 系 


也 们 的 模型 成 功 地 再 现 了 银 


河 系 中 心 核 球 的 垂 向 金属 丰 度 梯度 。 对 此 ， 他 们 的 解释 是 ， 由 于 银河 系 动力 学 演化 中 的 雅 死 


比 能 量 (Bj = E — Q L) 近似 守恒 ， 尽 管 存 在 强烈 的 动力 学 混合 过 程 ， 但 这 些 不 稳定 过 程 


在 混合 过 程 中 仍然 优先 占据 核 球 外 


属 恒星 主导 ， 从 而 在 观测 上 显示 出 垂 


并 不 会 改变 贫 富 金属 恒星 的 雅克 比 能 量 的 相对 分 布 。 这 会 导致 初始 位 于 外 


围 的 区 域 ， 


而 靠近 盘面 的 部 分 仍然 由 初始 位 


围 的 贫 金 属 恒星 


于 中 心 富 金 


名 的 金属 丰 度 梯度 。 这 种 解释 并 没有 得 到 广泛 的 认同 ， 
其 中 存在 的 问题 之 一 是 ， 在 银河 系 核 球形 成 的 过 程 中 ， 随 棒 共 转 的 引力 势 场 (02,) 不 能 保持 
较 稳 定 的 状态 ， 也 就 不 满足 雅克 比 能 量 近似 守恒 的 条 件 。 因 此 ， 我 们 认为 这 种 观点 不 足以 解 
释 银河 系 核 球 垂 向 金属 丰 度 梯度 形成 中 最 根本 的 机 制 。 

对 银河 系 核 球 中 垂 向 金属 丰 度 梯度 形成 的 另 一 种 解释 ， 


由 Bekki 和 Tsujimoto 呈 最 先 提 


出 ， 即 银河 系 核 球 并 非 由 单个 星系 盘 形成 ， 而 是 由 多 个 具有 不 同 金属 丰 度 分 布 的 星系 盘 ( 注 


盘 、 厚 盘 ) 共同 经 过 长 期 缓慢 演化 而 来 。 由 于 薄 盘 、 厚 盘 的 天 
星 平均 金属 丰 度 较 低 且 较 为 年 老 ， 垂 向 速度 弥散 更 大 ， 注 盘 的 恒星 更 年 轻 但 平均 金属 丰 度 较 


A] » 重 向 速度 弥散 较 小 ， 两 者 £ [n] 


BAL SEAN, EHSERS RE 


演化 形成 的 核 球 结构 中 厚 盘 恒 星相 比 薄 盘 恒 星 整 体 上 占 


据 更 大 的 垂直 空间 范围 ， 从 而 产生 垂 向 的 金属 丰 度 梯度 。Fragkoudi 4 A" 及 Di Matteo 等 


3 期 


AP 


种 观点 很 重要 的 一 个 观测 证 据 是 银河 系 核 球 的 金 
个 峰值 ， 这 意味 着 在 核 球 中 可 能 存在 多 个 星 族 的 共 


TERUS E. 


Debattista 等 人 "分 析 了 一 种 不 同 的 多 盘 模 型 
生 不 同 金属 丰 度 不 同 空间 分 层 的 主要 因 
[的 星系 盘 模 型 ， 且 初始 径 向 速度 弥散 随 着 金属 丰 度 的 增 


刘 仲 ， 等 ,银河系 核 球 中 垂 向 金属 丰 度 梯度 的 一 种 可 能 的 动力 学 起 源 机 和 


素 ， 并 假设 不 同 盘 的 初始 垂 向 速度 弥散 是 相同 的 。 在 
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= 


也 利用 类 似 的 具有 多 稳 结 构 的 多 体 模 拟 模型 来 解释 核 球 的 垂 向 金属 丰 度 梯度 的 形成 。 这 


二 


届 丰 度 分 布 直方 图 在 贫 富 金属 
司 演化 ， 而 这 些 星 族 成 分 极 有 可 能 来 自 


区 域 存在 多 
薄 


4， 他 们 认为 初始 的 径 向 速度 弥散 才 是 产 


他 们 的 多 体 模 拟 中 ， 使 用 了 五 个 于 力 
加 而 减 小 (同一 星系 盘 的 所 有 恒星 有 固定 的 金属 丰 度 )。 这 些 


Ei 
X 果 是 ， 


aN 


星系 盘 共同 演化 形成 的 结 


初始 径 向 速度 弥散 更 大 的 贫 金 属 盘 在 演化 中 呈现 出 更 大 的 垂直 空间 范围 ， 而 初始 径 向 
弥散 较 小 的 富 金属 盘 更 倾向 于 占据 盘面 附近 的 空间 。 这 表明 除 垂 向 速度 弥散 外 ， 星 系 盘 
向 速度 弥散 在 金属 丰 度 分 布 演化 中 
径 向 金属 丰 度 梯度 如 何在 长 期 缓慢 演化 模型 中 形成 垂 向 金属 丰 度 梯度 ”以 及 径 向 和 垂 
向 速度 弥散 在 金属 丰 度 演化 中 起 到 什么 样 的 作 
基础 上 提出 了 一 种 关于 银河 系 核 球 中 垂 向 金 


pu 


们 首先 利用 一 个 简 
BE (第 加 章 )， 并 描绘 出 模型 中 金 


我 们 通过 化 学 - 动力 学 的 分 析 来 说 明 核 球 中 垂 向 金属 丰 度 梯度 的 形成 过 程 ， 
何 形成 这 种 较为 


化 图 景 中 的 不 稳定 性 过 程 是 如 


可 能 占据 重要 地 位 ， 甚 至 是 主导 作用 。 


速度 
的 径 


]? 围绕 这 两 个 关键 问题 ， 我 们 在 前 人 研究 的 

属 丰 度 梯度 的 动力 学 起 源 更 为 合理 的 解释 。 我 

单 的 银河 系 核 球 的 化 学 -动力 学 模型 来 重 现 银河 系 核 球 的 垂 向 金属 丰 度 梯 
属 丰 度 在 核 球 中 的 整体 分 布 情况 (FB 章 )。 以 此 为 基础 ， 
解释 长 期 缓慢 演 


序 的 特征 的 (第 团 章 )。 另 外 ， 我 们 还 将 讨论 


这 一 模型 所 包含 的 创新 意义 和 可 能 存在 的 不 足 (第 回 章 )。 


2 化 学 -动力 学 模型 


本 工作 使 用 Shen 等 人 ”建立 的 多 体 数值 模拟 模型 。 在 这 个 银河 系 核 球 的 动力 学 模型 中 ， 


包含 1 x 106 个 模拟 粒子 的 
学 长 期 缓慢 演化 过 程 ， 即 星系 棒 


星系 盘 (标高 标 长 分 别 为 0.2 kpe 和 1.9 kpc) 通过 非 


AL, Ate yb 


常 简洁 的 动力 


的 不 稳定 性 


成 和 核 球 的 屈曲 不 稳定 性 过 程 ， 产 生 了 能 


很 好 地 匹配 BRAVA 巡天 ”中 关于 银河 系 核 球 动力 学 分 布 的 观测 结果 。 这 个 模型 也 能 很 好 地 


将 银河 系 核 球 中 的 


他 细节 结构 模拟 出 来 ” 。 我 们 以 Shen 等 人 的 数值 模拟 模型 为 基础 ， 并 


参照 了 Martinez-Valpuesta 和 Gerhard 等 人 呈 呈 使 用 的 化 学 -动力 学 建 模 方法 进行 了 研究 ， 


赋予 初始 状态 下 的 
ANA: 


星系 盘 粒 子 金属 丰 度 分 布 


， 设 径 向 金属 丰 度 作为 初始 的 金属 丰 度 梯度 ， 表 
[M/H]n = |M/H]o t on X R , (1) 


= 


KF, [M/H] 为 模拟 粒子 金 居 


E 
Py 


径 向 金属 丰 度 梯度 所 容许 的 范围 
始 的 垂 向 金属 丰 度 梯度 j 


=, [M/H] 为 模 
向 金属 丰 度 梯 度 (0.3) R 的 单位 为 kpc。 这 里 使 用 的 初始 
“由 内 而 外 ”(inside-out formation) 形成 


图 景 的 
(—0.15 ~ 一 0.3)”。 此 外 ， 考 虑 到 初始 
` 会 显著 影响 演化 结 


型 中 心 处 的 金属 丰 度 分 布 (0.6). ar AE 
径 向 金属 丰 度 梯度 是 星系 演化 
日 在 高 红 移 星 系 中 观测 到 的 星系 
星系 盘 非 常 薄 ， 初 
观测 和 理论 上 也 没有 初始 薄 盘 存在 高 


自然 结果 ， 
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垂 向 金属 丰 度 梯度 的 证 据 。 
为 0。 


因此 我 们 对 模型 做 了 进一步 简化 ， 将 初始 垂 向 金属 丰 度 梯度 设 


由 于 设 模 拟 粒 子 在 演化 过 程 中 的 金属 丰 度 值 不 变 ， 我 们 的 化 学 -动力 学 模型 中 的 化 学 分 


布 的 演化 完全 来 自 于 动力 学 过 程 带 来 的 模拟 粒子 的 重新 分 布 。 因 此 ， 这 个 模型 能 直接 地 回答 


在 模拟 星系 的 长 期 缓慢 演化 的 


不 稳定 性 过 程 中 不 同 金属 丰 度 的 粒子 是 如 何 进 行 混合 的 ， 以 


及 核 球 的 垂 向 金属 丰 度 梯度 能 否 在 这 个 过 程 中 形成 的 问题 。 


3 结 R 


3.1 ”金属 丰 度 分 布 直方 图 


为 了 去 除 星 系 盘 前 景 和 背景 的 盘 星 粒 


在 Ra < 4.5 kpe 范围 内 的 粒子 。 图 


子 对 于 核 球 观测 的 和 干扰， 我 们 选用 距离 银 心 半径 


四 展示 了 化 学 -动力 学 模型 中 在 核 球 范围 内 不 同 区 域 的 金 


属 丰 度 分 布 直 方 图 。 在 模型 中 核 球 区 的 不 同 银 纬 = —49. b= -—6°. b — 8 及 b= 二 —10? 


处 ， 我 们 用 黑色 箭头 标识 出 该 


区 域 金属 丰 度 分 布 的 平均 值 。 图 1 a) 为 模型 演化 时 间 2.4 Ga 


的 结果 ， 对 应 于 核 球 内 部 动力 学 演化 已 接近 完成 的 时 刻 ， 而 图 1 b) 为 演化 时 间 4.8 Ga 的 结 


果 ， 对 应 于 星系 盘 在 更 长 期 的 稳定 演化 的 结果 。 结 果 显 示 在 —8° <b < —4 范围 内 ， 核 球 
沿 着 短 轴 方向 的 垂 向 金属 丰 度 梯度 为 “0.039 dex/(e)， 与 Gonzalez 等 人 四 在 相应 区 域 的 观 
测 值 -0.04 dex/(°) 和 Gaia-ESO™ 的 观测 值 -0.05 dex/(°) 符合 较 好 。 另 外 ， 不 同 演化 时 


这 表明 银河 系 的 长 期 动力 学 演化 中 形成 垂 
中 的 不 稳定 性 演化 ， 即 棒 形 成 的 不 稳定 性 过 程 和 核 球形 成 的 屈曲 不 稳定 性 演化 。 
直方 图 中 显示 出 的 另 一 个 特征 是 其 分 布 的 非 正 态 性 。 前 文 提 到 银河 系 核 球 被 认为 可 能 


刻 金属 丰 度 梯度 分 布 的 对 比 也 显示 出 模型 所 产生 的 银河 系 核 球 中 的 金属 丰 度 梯度 是 稳定 的 ， 


名 金属 丰 度 梯度 的 主导 过 程 ， 是 核 球 在 形成 过 程 


具有 多 盘 复 合演 化 的 形成 历史 ， 其 重要 证 据 是 金属 分 布 直方 图 在 不 同 金属 丰 度 处 存在 多 个 
峰值 。 虽 然 贫 富 金属 丰 度 处 所 呈现 出 的 峰值 比例 在 观测 上 仍然 存在 争议 一 ， 但 是 均 呈 现 


出 非 正 态 的 直方 图 分 布 且 整 体 向 富 金属 端 偏 移 。 观 测 中 的 这 种 现象 通常 被 认为 是 银河 系 核 


球 具 有 复杂 结构 ( 薄 盘 、 厚 盘 、 尝 ) 共同 演化 的 证 据 ”。 对 此 ， 我 们 采用 了 IDL 的 xgaussfit 


程序 ”对 模型 结果 的 金属 


F 度 分 布 直 方 图 


解 ， 并 将 拟 合 的 高 斯 曲线 绘制 在 图 四 中 。 
态 特征 的 金属 丰 度 直方 图 ， 而 且 产生 的 贫 语 金属 峰 的 相对 比例 也 与 最 新 的 Blanco DECam 


Bulge Survey (BDBSS ) 巡天 观测 


进行 高 斯 混合 模型 (Gaussian mixture model) 分 


我 们 单 盘 的 化 学 -动力 学 模型 不 仅 能 够 产生 非 正 


的 金属 丰 度 分 布 结果 相符 合 。 对 于 银河 系 核 球 中 这 部 


分 贫 金 属 峰 的 起 源 ， 我 们 将 在 第 下 2 节 中 进行 讨论 。 由 于 不 同 观 测 中 金属 分 布 直方 图 的 结果 


目前 关于 银河 


系 核 球 中 是 否 存在 复杂 结构 共同 演化 仍然 是 一 个 


开放 性 问题 ， 厚 盘 和 党 的 部 分 在 核 球 中 的 分 布 可 能 并 不 会 非常 明显 。 


3.2 (1,5)- 金 属 丰 度 分 布 图 


除了 通过 金属 丰 度 分 布 直方 图 
和 银 纬 (5) 坐标 轴 上 的 金属 丰 度 分 布 ( 见 图 


外 ， 我 们 还 描绘 了 模型 中 所 模拟 的 银河 系 核 球 在 银 经 (D) 


p). 


星 比例 


日 


= 


Xip, S: RA RARER Pe I) oe JE BE — FT BB BL 


0.20 mi mi 
| (0,4) (0,74) 
"C cfl. | 
0.20 n^ E mmm 
(0,-6) (0.76) 
a —— emm. | 
(0,-8) (0,-8) 
0.10 
0.20 
(0,-10) (0,-10) 
0.10 
nae 15 -1.0 -05 00 05 10 -15 -1.0 -05 00 05 10 
[M/H M/H] 
a) b) 
FE: a) 模型 在 演化 时 间 为 2.4 Ga 时 的 分 布 直方 图 ，b) 模型 演化 时 间 为 4.8 Ga 时 的 分 布 直 方 图 。 黑 色 箭 头 标 
识 此 区 域 的 平均 金属 丰 度 值 ， 而 蓝 色 和 绿色 箭头 分 表 标 识 Gaia-ESO 和 Zoccali 等 人 图 在 此 区 域 观测 结果 的 


平均 值 。 高 斯 曲线 是 采用 高 斯 混合 方法 (GMM) 对 模拟 得 到 的 直方 图 进行 拟 合 的 结果 。 


1 “化 学 -动力 学 模型 对 银河 系 核 球 金属 丰 度 分 布 的 模拟 结果 


-0.06 
-0.12 
-0.18 
-0.24 
-0.30 
-0.36 
-0.42 
-0.48 
- -0.54 
) 


注 ， 两 条 虚线 中 间 代 表 银河 系 核 球 中 容易 被 尘埃 消光 而 难以 观测 的 部 分 。 
2 ”银河 系 核 球 模型 模拟 的 银 心 方向 (1, b) 轴 上 的 金属 丰 度 分 布 图 


45 b/(°) 


金属 丰 度 


5 0 
银 经 1/(° 
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为 了 与 观测 相 比 较 ， 我 们 选取 银 经 在 —10° < 1 < 10° HRH 
范围 内 ， 与 太阳 距离 为 4 ~ 12 kpe AMR. Enn 
距 银 心 8.5 kpc 处 (参考 文献 回 中 的 
中 几乎 在 各 个 方向 上 都 存在 一 定 的 金 忆 
而 核 球 的 中 心 部 分 仍然 由 富 金 属 丰 度 
的 趋势 相 一 致 。 进 一 步 地 ， 我 们 还 发 现 (1,5)- 金 属 丰 度 分 布 
现象 ， 体 现 为 1= 5° 处 的 金属 
测 结果 中 也 有 体现 。 对 此 ， 比 较 合理 


F l= 一 5° 的 相应 
的 解释 是 银河 系 棒 的 1 


41 卷 


在 —10° < b < 10° 核 球 
置 设 定 于 与 棒 主 轴 夹 角 20°， 
图 1)。 从 结果 中 可 以 看 到 ， 在 我 们 的 化 学 -动力 学 模型 
丰 度 梯度 ， 且 核 球 的 外 围 由 
的 模拟 粒子 所 主导 。 这 个 结果 与 VVV 观测 一 中 得 到 
图 中 存在 金属 丰 度 分 布 非 对 称 的 
区 域 ， 这 在 Gonzalez 等 人 "的 观 
硕 角 导 致 视觉 效应 ，1 > 0° 一 端 离 


盟 的 模拟 粒子 占据 ， 


观测 者 更 近 ， 而 1 < 0° 一 端 离 观测 者 较 远 ， 使 得 其 在 银 经 方向 展开 的 角度 不 同 。 这 些 与 观 


测 结果 相符 的 特征 ， 说 明 银 河 系 棒 和 
属 丰 度 分 布 的 演化 起 到 非常 重要 的 作用 。 


3.3 ”化 学 -动力 学 关系 


银河 系 核 球 中 不 同 金属 


核 球 的 化 学 -动力 学 演化 对 于 研究 和 


lLE 度 的 星 族 动力 学 特性 


核 球 中 经 历 星系 盘 长 期 缓慢 演化 的 
机 速度 弥散 ， 而 由 其 他 动力 学 演化 过 下 


族 应 该 保留 更 多 的 盘 状 特 行 


区 域 的 金 


成 演化 历史 的 重要 因素 。 
E， 如 柱状 旋转 和 较 低 的 随 
BEER) 形成 的 星 族 则 会 呈现 与 之 


不 同 的 动力 学 和 化 学 特征 ， 如 较 高 的 随机 速度 弥散 。Ness 等 人 考虑 根据 金属 丰 度 的 不 同 
将 核 球 中 的 恒星 分 解 为 [Fe/H] > 0 (A), —0.5 < [Fe/H] <0 (B) 41 —1.0 < [Fe/H] < —0.5 


(C) 三 个 主要 星 族 和 [Fe/H] < —1.0 B D 星 族 。 他 们 发 现 三 个 主要 星 族 : 
柱状 旋转 ， 但 在 速度 弥散 的 分 布 上 富 金 属 的 


状 特征 的 
族 A 和 B 表现 出 较 强 的 垂 向 梯度 ， 而 贫 金 属 


的 星 族 C 在 垂直 方向 的 变化 中 ， 其 速度 弥散 只 产生 很 小 的 差异 。 而 我 们 的 模拟 中 ， 不 同 金 


属 丰 度 的 星 族 除了 随机 速度 弥散 随 着 金属 
速度 弥散 的 垂 向 分 布 梯度 ) 并 没有 显著 的 


FE 度 降低 而 升 高 外 ， 动 力学 基本 特 和 


E (柱状 旋转 、 


T A) [Fe/H]>0 
a o (A) [Fe/H] 
B 
c 0 
Phe 
o 
© -50 
-100 
120 
p^ 
4 80 
ao -F oa 
5 40} 
-10 -5 


: 带 误 差 棒 的 方块 和 点 代表 ARGOS 中 的 


型 中 在 相应 区 域 的 模拟 结果 (黑色 为 5 — —5, HEJ b= 一 10°)。 


图 3 ”银河 系 核 球 化 学 -动力 学 模型 的 不 同 星 族 的 旋转 曲线 (上 ) 和 速 


10*)， 实 线 代表 模 


3 期 刘 仲 ， 等 ; 银河 系 核 球 中 垂 向 金属 丰 度 梯度 的 一 种 可 能 的 动力 学 起 源 机 于 


如 果 贫 金属 丰 度 星 族 的 随机 速度 弥散 分 布 的 差异 确实 显著 存在 ，Ness 等 人 ”提出 ， 


金属 丰 度 在 -1.0 < [Fe/H] < -0.5 范围 内 ， 银 河 


化 历史 的 厚 盘 结 构 。 但 是 ， 最 近 Wylie 等 人 采用 Ness 等 人 的 ARGOS Side", FFA 
APOGEE 的 观测 结果 ”” 进行 了 交叉 匹配 ， 发 现 星 族 C (—1.0 < [Fe/H] < —0.5) 的 速度 弥 
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= 


系 核 球 更 有 可 能 属于 与 薄 盘 具有 不 同 演 


散 分 布 没 有 星 族 A 和 星 族 B 的 明显 。 我 们 认为 Ness 等 人 “所 认定 的 星 族 C JT I dile 
的 观点 相应 支持 的 证 据 并 不 十 分 充分 ， 且 银河 系 核 球 中 不 同 金属 丰 度 的 星 族 可 能 主要 仍然 


是 由 薄 盘 长 期 缓慢 演化 而 来 。 


4.1 ， 核 球 垂 向 金属 丰 度 梯度 的 形成 机 制 


度 的 形成 是 源 于 核 球形 成 过 程 中 的 两 种 不 稳定 性 ， 


根据 之 前 的 分 析 和 讨论 ， 我 们 知道 银河 系 核 球 的 化 学 -动力 学 模型 中 的 垂 向 金属 丰 度 梯 


径 向 和 垂 向 的 速度 弥散 在 这 个 过 程 中 


起 到 非常 重要 的 作用 。 我 们 分 析 不 同 金属 丰 度 的 星 族 在 模型 不 同 演化 时 刻 的 径 向 和 垂 向 速 
度 弥 散 分 布 ， 结 果 如 图 4 所 示 。 它 们 分 别 是 初始 状态 的 了 = 0 时 刻 ， 星 系 棒 结构 形成 但 未 
开始 届 曲 不 稳定 演化 的 了 = 0.6 Ga 时 刻 ， 形 成 核 球 的 屈曲 不 稳定 性 演化 过 程 中 (了 = 1.8 


Ga 时 刻 ， 接 近 演 化 完成 )， 核 球 结构 演化 完成 并 达到 稳定 后 (T= 4.8 Ga 时 刻 )。 不 同 于 图 


男 ， 本 章 我 们 采用 了 更 加 细致 的 星 族 分 层 ， 根 据 金属 丰 度 的 不 同 将 参与 核 球形 成 的 恒星 分 为 
0.3 ~ 0.6 dex (初始 半径 : 0~ 1 kpc), 0 ~ 0.3 dex (初始 半径 : 1 ~ 2 kpc), —0.3 ~ 0 dex (XJ 
始 半 径 : 2~ 3 kpe), —0.6 —0.3 dex (初始 半径 : 3 ~ 4 kpc) 和 一 0.9 ~ —0.6 dex (初始 
半径 ，4 ~ 5 kpc) 五 个 成 分 ， 并 绘制 出 其 径 向 和 垂 向 速度 弥散 的 演化 。 


图 团 a) 显示 了 模型 在 初始 状态 下 不 同 星 族 的 径 向 速度 弥散 or 和 垂 向 速度 弥散 oz (ER 


= 5 kpe 内 的 分 布 ， 分 别 用 不 同 颜色 的 实 线段 和 虚线 段 表 示 。 可 以 清晰 地 看 到 不 同 半径 处 的 


金属 丰 度 和 相应 的 平均 速度 弥散 ， 随 着 半径 的 增加 而 降低 。 这 是 典型 的 星系 盘 密 度 指 数 分 布 


模型 的 速度 弥散 分 布 特征 。 图 田 b) 显示 星系 棒 的 形成 造成 了 模拟 粒子 强烈 的 径 向 混合 ， 不 
同 星 族 在 垂 向 速度 弥散 oz 几乎 不 变 的 同时 ， 径 向 速度 弥散 op 产生 显著 的 变化 。 此 时 贫 


属 星 族 特别 是 进入 盘 中 心 2 kpe 范围 的 部 分 其 径 向 


速度 弥散 or 急剧 增加 ， 与 之 相 比 富 金属 


星 族 的 径 向 速度 弥散 or 几乎 没有 发 生 明 显 的 变化 。 当 发 生 核 球形 成 的 屈曲 不 稳定 性 演化 时 


( 见 图 团 c))， 不 同 星 族 的 径 向 速度 弥散 or 并 不 发 生 明 显 的 演化 ， 而 垂 向 速度 弥散 oz 显著 
增加 ， 且 这 个 过 程 主要 发 生 在 中 心 的 核 球 区 域 (0 ~ 2 kpc)。 核 球 届 曲 不 稳定 性 演化 的 最 终 


结果 是 ， 垂 向 速度 弥散 和 径 向 速度 弥散 的 比值 oz/on 会 达到 一 个 较为 稳定 的 值 ”。 因 此 ， 
在 棒状 不 稳定 性 演化 中 获得 了 较 高 的 径 向 速度 弥散 op 的 贫 金 属 星 族 ， 在 届 曲 不 稳定 性 演 


化 中 会 获得 较 高 的 垂 向 速度 弥散 oz， 且 这 意味 着 这 些 星 族 在 空间 上 会 占据 更 大 的 垂直 范围 。 
我 们 测定 的 垂 向 速度 弥散 与 径 向 速度 弥散 的 比值 cz/ca 约 为 0.8， 这 远 高 于 理想 薄 盘 的 理 
论 值 oz/or ~ 0.3， 但 在 更 接近 于 实际 的 星系 数值 模拟 中 是 比较 合理 的 ( 见 参考 文献 2 )。 


并 且 在 模型 中 的 核 球形 成 之 后 的 演化 过 程 中 ， 这 村 


的 垂 向 速度 弥散 与 径 向 速度 弥散 的 比值 
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一 0.3<[M/H]<0.6 
= 0«[M/H]«0.3 
200} | = -0.3<[M/H]<0 
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a) b) c) d) 


VE: 实 线 代 表 径 向 速度 弥散 wm， 虚线 代表 垂 向 速度 弥散 oz 的 分 布 情况 。 从 a)、b)、c)、d) 分 别 对 应 于 模型 演 


化 的 不 同 阶段 ， 分 别 为 模型 初始 状态 T = 0， 棒 结构 形成 初期 的 全 = 0.6 Ga， 核 球 屈曲 不 稳定 性 演化 过 程 中 
T = 1.8 Ga， 及 核 球形 成 后 经 过 长 期 缓 变 演 化 达到 全 = 4.8 Ga. 


图 4 


银河 系 核 球 的 化 学 -动力 学 模型 的 不 同 星 族 在 不 同时 刻 的 动力 学 参数 


cz/og 也 将 几乎 保持 稳定 ( 见 图 男 d))， 这 代表 核 球 的 垂 向 金属 丰 度 梯度 也 会 保持 稳定 。 


我 们 把 银河 系 核 球 中 金属 丰 度 分 布 的 形成 机 制 总 结 为 “两 步 ” 过 程 。 
的 不 稳定 形成 过 程 使 部 分 初始 半径 较 大 ， 金 属 丰 度 较 低 的 恒星 进入 中 心 
区 域 ， 并 且 这 部 分 恒星 获得 了 较 大 的 径 向 速度 弥散 ono 


1) 银河 系 中 棒 


2) 核 球形 成 的 


LE 


"ek. PEARLS Y BOKATSISIBEBGESE URBC or 的 贫 金 属 恒 


星 再 次 获得 较 大 的 垂 向 速度 弥散 wz， 从 而 占据 更 大 的 垂直 空间 范围 ， 而 核 球 中 径 向 速度 弥 


散 og 较 小 的 富 金 忆 


属 丰 度 梯度 。 


因此 ， 在 银河 系 核 球 的 垂 向 金属 


粒子 仍然 留 在 中 心 靠近 盘面 区 域 。 在 这 一 过 程 中 核 球 内 形成 了 垂 向 的 金 


E 度 梯度 的 形成 中 ，Debattista 等 全 讨论 的 径 向 速 


x 


度 弥散 与 Bekki 和 Tsujimoto"" 所 讨论 的 垂 向 速度 弥散 均 在 演化 中 起 到 重要 作用 ， 并 且 两 者 


通过 核 球 演 化 形成 的 相对 固定 的 比值 相 联系 。 核 球 中 垂 向 速度 弥散 与 径 向 速度 弥散 的 比值 
az/aaR 说 明星 系 盘 的 内 部 演化 在 影响 各 星 族 径 向 速度 弥散 分 布 的 同时 ， 也 会 影响 其 最 终 的 


径 向 速度 弥散 的 分 布 。 我 们 认为 ， 这 种 化 学 -动力 学 演化 机 制 应 当 也 存在 于 其 他 的 银河 系 长 
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Abstract: Observations show a clear vertical metallicity gradient in the Galactic bulge. It 
is often regarded as a counter-argument against the ‘secular evolution’ of the Milky Way 
bulge, as the strong mixing of the internal instabilities in the secular evolution may erase 
the pre-existing vertical chemical gradient. Inspired by some previous work, we use a thin 
disk model that self-consistently forms a boxy/peanut bulge via the bar and buckling insta- 
bilities, to explore the origin of this chemical gradient. We find the vertical gradient can be 
established if the initial disk is set up with a radial metallicity gradient, as in the evolution 
the metal-poor particles in the outer region move inward and get more vertically thickened 
than the inner metal-rich particles. This is because the formation of the bar always tend 
to heat the in-plane random motions (og) of the metal-poor population at larger radius 
more, therefore when the oz/or of all populations finally reach a similar value about 0.8 
in the buckling instabilities, the vertical random motion (oz) of metal-poor population is 
also more heated. This formation mechanism revealed by our simplest secular evolution disk 
model must also be in action in other more complicated disk models that experienced similar 
internal instabilities. Our chemo-dynamical model can produce the multiple ^peaks" in the 
metallicity distribution functions found in many recent bulge observations. We also discuss 


other implications and the weaknesses of this model. 


Key words: Galaxy; Galaxy bulge; Galaxy evolution; Galaxy abundance 


